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Abstract 

We have studied the integer quantum Hall effect at finite temperatures 
by diagonalizing a single body tight binding model Hamiltonian including 
Aharonov-Bohm phase. We have studied the energy dependence of the spe- 
cific heat and the Hall conductivity at a given temperature. The specific heat 
shows a sharp peak between two consecutive Hall plateaus. At very low tem- 
peratures, the numerical results of the temperature variations of specific heat 
(in the plateau region) are in good agreement with the analytical results. 
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Quantum Hall effect continues to draw considerable interest since its surprised discovery 
|jlj] in two dimensional electronic systems at very low temperatures and in the presence of 
very high magnetic fields. In this phenomenon, the Hall conductivity gets quantized at 
integer filling factors v with the value ve 2 /h in high accuracy and the diagonal conductivity 
vanishes at the same time. A similar phenomenon also occurs at certain fractional filling 
factors as well f2|. In this paper, however, we are concerned only with the integer quantum 
Hall effect (IQHE). 

The scaling theory [|J predicts that all the states are localized in two dimensions due to 
disorder potential which respects time reversal symmetry. In the presence of strong magnetic 
field, perpendicular to the plane of the system, the electronic energy spectrum of the system 
shows a series of Landau bands. Since the magnetic field destroys time reversal symmetry 
and suppresses back scattering, some of the states appear to be extended at the center of 
each Landau band. On the other hand, the states away from the center of each band remain 
localized. The transport properties (electronic conduction) of the system do depend on the 
distribution of localized and extended states. This distribution of localized and extended 
states gives rise to IQHE. So long as the Fermi level lies in the region of localized states, the 
Hall conductivity remains pinned to a value forming a plateau. The critical transition || 
from one plateau to another occurs when Fermi level passes through the region of extended 
states. 

On the other hand, the equilibrium properties (say specific heat) of the system depend 
simply on the density of states These properties cannot distinguish the localized 

states from extended states. This is indeed so, as it emerges from our systematic numerical 
diagonalization study in a lattice model. Agreeably, the separation of electronic specific 
heat from lattice specific heat in quantum Hall systems, from the experimental data, are 
extremely difficult. However, the experiments 0,^] are performed in multilayer systems to 
extract electronic specific heat. 

In this paper, we consider a single body (spinless) tight binding model Hamiltonian || 
with random impurities in each site in the presence of a uniform magnetic field perpendicular 
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to the lattice plane. We have calculated the Hall conductivity a H and the specific heat C v 
at a finite temperature T in a finite lattice. We show the formation of Hall plateau and the 
transition between two plateaus of the Hall conductivity by studying the energy variation 
of the Hall conductivity. When the Fermi level lies between two Landau bands, C v is 
exponentially small since in such a case the probability of low energy excitation is very 
small. In such a region we find C v ~ (l/T 2 )e~ A//T , where T is the temperature and A is 
a typical interband excitation gap. This behavior of C v is consistent with the analytical 
results obtained from the self consistent Born approximation. On the other hand, when 
Fermi level lies within a Landau band C v increases substantially and shows a sharp peak. 
This peak position of C v lies between two consecutive Hall plateaus. We show analytically 
that C v in this region is, in fact, proportional to the density of states at Fermi energy. C v 
reflects the profile of the density of states. The width of the Hall conductivity plateau is 
larger than that of C v plateau. This is due to the fact that while the transport properties 
distinguish between localized and extended states, the equilibrium properties do not. 
Single body two dimensional tight binding Hamiltonian is given by 

H = ~ E (' ( ') ( \, + h.c) + E (!) 

<ij> i 

where Cj and Cj are spinless Fermionic creation and destruction operators respectively at 
lattice site i. Here < ij > refers to nearest neighbour hopping, and the hopping parameter 
is taken to be unity. A uniform magnetic field B is applied such that it produces a uniform 
magnetic flux = (hc/e)/M (M integer) per plaquette. When a particle hops from site j 
to site i, it acquires a phase = J] A ■ dl, where the vector potential in the Landau gauge 
is taken to be A = (0, Bx). For this choice of gauge, an electron gets nonzero phase while 
hopping along y direction only. The (diagonal) second term in the Hamiltonian corresponds 
to random disorder potential at site i with strength Wi varying from — W/2 to W/2. We 
have constructed the Hamiltonian (1) by considering the periodic boundary condition. As 
a result, the energy spectrum becomes M number of Landau bands. Each band contains 
L 2 /M states in a square lattice of size L x L. 
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We compute the Hall conductivity ou using the Kubo formula |J 

ze 2 ft x ^ x ^ (/ n - f m ) < m\v x \n >< n\v y \m > 

where E n is the energy of |n >-th state, and 



i/h) ]T [ ( ']-: ( 'i< • J - tfCj+re-***™ ) (3) 



is the velocity operator fl(J along direction r = (x or And f n = l/(exp[(E n —Ep)/kBT] + 



1) is the Fermi function at temperature T with Fermi energy Ep. The electronic specific 



heat per unit area of the system can be computed using the formula [11 



C v = — ^ £(£„ - E F ff n {\ - f n ) . (4) 

Figure 1 shows the energy variations of a# and C v . These are computed using the 
formulae (2-4) for a system of size 12 x 12, strength of random potential W = 0.5, flux 
per plaquette (ft = (hc/e)/12, and at small temperature k B T = 0.01. Three distinct IQHE 
plateaus, between four lowest Landau bands, at oh = u(e 2 /h) [y = 1,2,3) are shown 
for Ep < 0. It is observed that the higher plateaus are formed at relatively lower values 
|T2| | compared to corresponding integral multiple of e 2 /h. The decrease in the value of 



quantization is due to the relatively higher values of impurity strength for the formation of 
the corresponding plateaus Jl2| . The transitions between plateaus occur at the centers of the 
Landau bands. In these transition regions, C v shows sharp but finite peaks with the width 
same as the Landau band width. In between the Landau bands, C v is exponentially small 
since the low energy excitations are the least probable and the plateaus in C v are also formed. 
The widths of plateaus in an are relatively higher than that of the corresponding plateaus 
in C v . This is because of the fact that the transport properties (like an) are sensitive to 
extended states or localized states which are formed at the tail of the Landau bands but the 
equilibrium properties (like C v ) are insensitive to the nature of the states. This is indeed 
the case since the expression (4) for the specific heat does not depend on the characteristics 
(localized or extended) of the states. 



Figure 2 shows the temperature variation of the specific heat at two different energies 
corresponding to two different plateaus (see Fig.l), keeping other parameters fixed. The 
best fitted form of C v is found to be C v ~ (l/T 2 )e~ A ' T which is consistent also with the 



numerical diagonalization result |TJ| for fractional quantum Hall states. Here A is a typical 
excitation gap which we find to be consistent with the Landau band gaps. The values of A 
for these two curves are found to be almost same while the prefactors are different. 

We now analytically evaluate the specific heat at the regions of the plateaus and the 
transition between the plateaus separately at low temperatures by using self consistent Born 
approximation (SCBA) for the impurity potential, and compare those with the numerical 



results. The density of states of an IQHE system in SCBA |14J is given by 

which is a series of semicircular bands of radius T at Landau levels with energy e n = (n + 
^)hu c . Here u c is the cyclotron frequency and / is the magnetic length of the system. We 
assume that the broadenning 2T of each Landau level is very small compared to %u c , i.e., 
T << Tuo c . The electronic specific heat of the system can be expressed as 

1 r°° 

C -=k^5j ~ E F) 2 fm - f(e)]de , (6) 

where /(e) = 1/ ( e [ ^' EF)/kBT] + 1) is the Fermi function. 

(i) Specific heat at the plateau regions: The specific heat can be expressed using eqns (5) 
and (6) as 

H r e n+ F r ll/2 

C^rWE/ {€ ~ Ef)2 r2 -( e " e ^ 2 f(e)[l-f(e)}de, (7) 

with the constant C = 1/(tc 2 1 2 T 2 ). We now consider that the Fermi energy Ep is in between 
two consecutive Landau bands. In this case an integer number of bands are filled. By 
Taylor's expansion of eqn (7) in powers of T /huj c , we obtain 
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Cv 



CitT 
'k B T 2 ~ 



J2(e n -E F ) 2 f(e n )[l-f(e n )} + 



r 



(8) 



At low temperatures (hu c /kBT » 1), we perturbatively evaluate C v with the expansion 
parameter JT3I] e -^c/ fe s T . We thus find 



C V (T) 



nl 2 k 



B 



u 2 N 2 



- E F u c N 



J^2 



-%uj c /2k B T 



(9) 



where N is the number of filled Landau bands. Note that the specific heat is activated 
through the gap between two Landau levels. Interestingly, this is the same temperature 
dependent form of C v that we have found from our numerical diagonalization study. The 
prefactor in Eqn (9) depends on the value of Fermi energy and the number of filled Landau 
bands. For this reason, we obtain two different curves for two different Fermi energies in 
Fig. 2. 



(ii) Specific heat at the regions of transition between two consecutive plateaus: In these 
regions, Fermi energy lies inside a Landau band. The approximations made in (z) will not 
hold good in this case since closeby states above the Fermi energy are available for low energy 
excitations. Using Sommerfeld's low temperature expansion |9| of the energy derivative of 
Fermi function in eqn (6), we obtain 

C v = ^k 2 B Tp(E F ) . (10) 

Clearly, the specific heat in the transition regions follows the profile of the density of states 
which has maxima at the centers of the Landau bands. This is indeed so that we have also 
obtained this type of profile from our numerical study. However, the linear temperature 
dependence that we have obtained in eqn (10) does not resemble with our numerical calcu- 
lations. In our numerical study, we have obtained a fitted form C v ~ (l/T 2 )e _A//T here as 
well, with A much less than that in the case when Fermi energy lies between two successive 
Landau bands. We remark that this fitted form, in our numerical study, is a consequence 
of finite size diagonalization study since in such a case, there always remains small but fi- 
nite gap between the states within a Landau band. In other words, the linear temperature 
dependence of C v in eqn (10), obtained from continuum description of the density of states 
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(Eqn (5)) within a Landau band, cannot be achieved by a finite size diagonalization study. 
However, we believe that in L — > oo limit, the linear temperature dependence of C v may 
hold good. 

In summary, we have studied the finite temperature Hall conductivity and the specific 
heat of an IQHE system by finite size diagonalization of a tight binding model Hamiltonian 
including Aharonov-Bohm phase. The energy variation of Hall conductivity and the specific 
heat have been studied at low temperatures. The specific heat also shows plateau together 
with the Hall conductivity, and gets sharply peaked between two consecutive plateaus. We 
have also calculated specific heat analytically in self consistent Born approximation. The 
temperature variation of specific heat (in the plateau region) agrees well with the numerical 
results. 
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Figure Captions 

Fig. 1: Energy variations of <j H (circles) and C v (triangles), for ksT = 0.01, W = 0.5 and 
M = 12. a H is shown in the unit e 2 /h while C v is in arbitrary unit. 

Fig. 2: Temperature variations of C v for two different energies corresponding to two different 
plateaus, (i) Circles represent Ep = —3.0 and (ii) Triangles represent Ep = —2.2. Solid lines 
represent the best fitted functional form C v ~ (l/T 2 )e~ A / T . Here W = 0.5 and M = 12. 
Here, C v is shown in arbitrary unit. 
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Fig. 1 
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